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ABSTRACE 


A ten-level primitive equation ocean circulation model 
is used to PR OMER: « the formation and evolution of 
large-scale thermal anomalies observed in the central North 
Pacific Ocean during the fall and winter of 1976. Several 
initial value model integrations of 4 months duration are 
carried out in order to help explain the observed anomaly 
development. Initial and verifying ocean data down to 400 m 
depth are obtained from the NORPAX ships of opportunity pro- 
Gram. Anomalous atmospheric wind forcing is obtained from 
Namias' monthly mean sea-level pressure anomalies, while 
Climatological heat fluxes are used. 

The skill with which the model simulates the observed 
anomaly evolution in the different experiments is estimated 
synoptically and measured statistically by calculating root 
mean square (RMS) temperature errors and S1 skill scores. 
Analysis indicates that anomalous atmospheric wind forcing 
improves the model predictions in the upper levels. For 
this particular winter case usSing climatological heating, 
however, knowledge of the initial anomalous temperature con- 
ditions does not improve the model results. The model skill 
at upper levels exceeds both persistence and climatology 
(forecast of zero anomaly) while at the lower levels it is 


comparable to persistence but not as good as climatology. 
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Bis eon se LGURES 


Fig. 1. Geographical area occupied by the idealized 
rectangular ocean model; the TRANSPAC (ocean) 
data area and the method of extrapolating 
the oceanic temperature data outward into 
the data void regions of the model along the 
line AO=-j-efefrrrrrererrrrrrrrereres 


Fig. 2. Anomalous pressure fields in the central 
Pacific for September 1976 to January 1977 
at 2 mb intervals; shaded areas represent 
negative anomalous pressure; area shown is 
the area east of the dateline in which the 
temperature data are given in figures 3 
ty REACT CPL De ee Sie re re ye ee eee ee 


Fig. 3. Observed temperature anomaly fields for 
September 1976, November 1976 and January 
lov) woatetOywelLoz, Loe and 262 meters; con- 
POUr Interval at LO meters 1s 0:;5°C; con= 
EOUr Interval at 102, 162 and 262 meters 
is 0.2°C; shaded areas represent negative 
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meog. 4. Same as figure 3 except for Model 1 results 
Fig. 5. Same as figure 3 except for Model 2 results 
Fig. 6. Same as figure 3 except for Model 3 results 
Fig. 7. Same as figure 3 except for Model 4 results 
Fig. 8. Same as figure 3 except for Model 5 results 
Fig. 9. Plots of root mean square error (RMSE) for 


persistence (P) and ocean models (1, 2, 3, 


4 and..5) nondimensionalized by the RMSE 
Eouwclimatology, as a function of time in 
Gaysmuror lO, 102, 162 and 262 meters =- = 
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I. INTRODUCTION 


Early investigations by Namias and others indicated that 
the winds are an important source of kinetic energy for the 
oceanic circulation. But just as the mean atmospheric cir- 
Culation helps drive the mean large scale ocean circulation, 
anomalous atmospheric forcing can be assumed to drive anoma- 
lous currents in the oc€ean (Malkus, 12962). This premise was, 
in part, the motivating force for the present model study. 
The effects of anomalous atmospheric forcing on the oceans 
have been a matter of intuitive guessing for years, but few 
Studies have attempted to assess in more than a descriptive 
Or qualitative manner the nature or magnitude of these 
effects. A notable exception is the recent work of Davis 
Wee76, 1977). 

One of the most important manifestations of the effects 
of anomalous atmospheric forcing on the oceans is the devia- 
tion of the ocean thermal structure from the long term cli- 
matological norm. The space and time scales of these devia- 
tlons or anomalies are known in avery general sense, but 
the exact generating mechanism or mechanisms remain unclear 
at present. heiesmele Intent Of this investigation, within 
the larger context of the Anomaly Dynamics Study (ADS), to 
"Describe and understand the mechanisms responsible for 
large scale climatological anomalies and related changes in 


eacethermal structure and the general circulation of the 





ipper layers OF the ocean in order to predict their forma- 
ereon and evolution.” (Kirwan, 1976) THe reenter Lout lon. of 
anomalous atmospheric wind forcing as a generating mechanism 
for changes in the thermal structure of the ocean will be 
emphasized. 

Previous model investigations of the dynamics of large 
scale temperature anomaly formation in the ocean have been 
limited by the lack of verifying data, at the surface and 
at depth, over sufficiently large areas. In cases where 
anomalous atmospheric forcing for a certain period of time 
maroeaddaed to climatological wind forcing (Namias, 1965; 
Arthur, 1966; Jacob, 1967), model results show an important 
influence on the ocean by the anomalous forcing. 

In the present study, an ocean circulation model is 
used to make a number of initial value predictions of the 
evolution of large scale ocean anomalies. The different 
model predictions are designed to examine the effect on the 
Mmmeairction of including or excluding certain phySical pro- 
cesses in the model. The study is unique in that a care- 
fully collected and analyzed set of ocean temperature data 
are used. 

It must be emphasized at this point that this study 15 
not meant to be a test of the predictive capability of the 
model itself. First, a very important part of the atmos- 
pheric forcing fields, namely anomalous heat flux that 1s 
required for the testing of the prediction capabilities of 


any ocean model is not available on the time or space scales 
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required. Consequently, only anomalous winds are used in the 
present studies. Second, ocean modeling itself is still 
relatively unsophisticated and certainly untested in model- 
ing large scale thermal anomalies observed in the oceans 

even if all the atmospheric forcing were known. In addzeron, 
this study will serve aS an initial investigation into the 
adequacy of currently available temperature data in describ- 
ing the true thermal condition of the Central Pacific. With 
the current oceanic and atmospheric data available and with 
the current state of the science of ocean modeling, all that 
can be tested are a number of hypotheses regarding the rela- 
tive importance of various dynamic forces that act to drive 
the oceans. Any conclusions resulting from this study then 
are subject to re-evaluation as ocean modeling and monitoring 
Capabilities improve. The 1976-77 fall-winter period was 
selected for this study because of the large amplitude anoma- 
lies in the atmosphere and ocean that existed in the North 


Paetric at that time. 
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a DATA PROCESSING 


The modeling experiments described below made use of 
anomalous atmospheric winds and oceanic temperatures. The 
fPerst step in preparing the atmospheric data for input into 
the ocean model was calculation of sea-level pressure anoma- 
lies or departures from normal for all of 1976 and January 
and February of 1977. All the atmospheric data were generous- 
ly made available to us by Dr. J. Namias of Scripps. Long 
term (20 year) mean values of sea level pressure for each of 
these months, as determined by Namias, were subtracted from 
monthly means in order to obtain the sea level pressure de- 
partures from norm. Analysis of these anomaly fields indi- 
cated the period from September 1976 until January 1977 ex- 
hibited sufficiently large anomalous sea level pressures to 
be suitable for our investigation. 

In order then to input Namias' anomalous pressure data 
into the model, it was necessary to extrapolate and then 
interpolate the values from Namias' 5° latitude by 5° longi- 
tude (9x19) fields encompassing 20°N to 60°N and 145°E to 
L25°W to the 33x33 fields extending from the equator to 65°N 
and 145°E to 125°W required by the model (see figure l). 

The extrapolation of the anomalous pressure data toward the 
northern boundary assumed a slope that decreased by 50 per- 
cent at each grid point. Extrapolation towards the equator 


was such that anomalous u and v components of the wind field, 
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calculated from the anomalous pressures, go to zero at the 
equator. Once the extrapolation was completed, data north 
of 50°N was smoothed several times in the east-west direc- 
tion in order to remove any small scale features introduced 
by the extrapolation. This procedure was carried out for 
the months of September 1976 through January 1977. 

With the 33x33 field of sea level pressure anomalies 
computed, the corresponding anomalous geostrophic u and _v 
components of the wind were calculated. These anomalous geo- 
strophic wind fields were the only anomalous atmospheric 
conditions used in the ocean model. These fields (figure 2) 
are described in the Results section below. 

The next step in data preparation for input into the 
model was to compute the temperature anomalies that would 
serve aS initial and verifying conditions for the model. 

All temperature data was provided by TRANSPAC, the ships of 
Opportunity cooperative research effort. The TRANSPAC ex- 
pendable bathythermograph (XBT) field program is a part of 
the North Pacific Experiment (NORPAX) which is charged with 
investigating how oceanic variability influences the weather 
and climate, and in turn how the weather and climate varia- 
Peebmty influence the oceanic circulations. Prior to TRANSPAC, 
data on the thermal structure of the ocean was inadequate, 
thus leading to the inception of the TRANSPAC XBT field pro- 
gram to provide the subsurface temperature data necessary 
for NORPAX. The observational phase of the program began in 


miowany L975 utilizing shaps of opportunity (i.e. merchant 
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Smmees Voluntarily providing a no-cost observation platform) 
that routinely traverse the sea lanes between the United 
States and Japan. The temperature anomaly data provided 
through TRANSPAC was obtained by subtracting eee 
Means from the monthly mean observed temperatures. Anomalous 
temperatures were provided at ten levels including 0,30, 60, 
meee 20, 150, 200, 250, 300 and 400 meters. For initializ- 
ing temperature fields in the model, the anomalous tempera-~ 
ture field for September 1976 was added to the model gener- 
ated climatology for September. Calculation of verification 
statistics entailed use of anomalous temperature fields for 
Setober through December of 1976 and January 1977. Again, 
the TRANSPAC anomaly data had to be processed for use in the 
model because the domain and standard levels of the TRANSPAC 
data- did not coincide with the domain and standard levels of 
the model. Tt was decided, in order to avoid major changes 
Eameeie Model, to extrapolate from the 2° latitude by 5° longi- 
tude (11x15) TRANSPAC fields extending from 30°N to 50°N 
GmaeLOO°E to 130°W to the area and gridpoints encompassed 

by the model (see figure l). 

The first step in preparing the TRANSPAC temperature anom- 
aly data was to interpolate the data from the ten levels of 
TRANSPAC to 11 intermediate levels of the model, specifically 
meee, 4550900, —b25, 200, 325, 600, 1200, 2200 and 4000 
meters. Below 400 meters, values were assigned to the inter- 
Mediate levels by decreasing the 400 meter anomalous tempera- 


ture to a value of zero at the 2200 meter intermediate level. 
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This was accomplished by multiplying the 400: meter value by 
a cosine function which decreased with depth. For example, 
if T' (400) is the given TRANSPAC anomaly at any grid point, 


maem the anomaly at depth Z, T’'(Z), was specified by 


22400. 1 


Ti 2) = £400) cos l (T3860? 3] : 


This is a very important assumption, since it assumes the 
temperature anomaly cannot change sign with depth below 
400 meters. The resulting temperatures at the 11 intermedi- 
ate levels were then averaged between any two levels to de- 
Fine the temperature at one ten prognostic levels in the 
Menem wikech are located at 10.0, 32.5, 62.5, 102.5, 162.5, 
meee, ¢€62.5, 900.0, 1700.0 and 3100.0 meters. 

The next step in preparing the TRANSPAC data was to ex- 
pand it to cover the entire area of the model. This was 
done by expanding the temperature anomaly at the border of 
the 11x15 TRANSPAC area outward to zero at the outer boundary 
Smmrne model. In order to provide a smooth transition into 
the area of no data, the extrapolation was along a sine func- 
tion as shown in figure 1. Areas of the array which could 
be filled by both longitudinal and latitudinal extrapolation 
were assigned a temperature anomaly value which was an 
average of the two extrapolations. Once the entire geographi- 
cal area of the model was filled on a 2° latitude by 5° 
longitude grid, an interpolation scheme was used to fill 
Bie mo oxaoearray OL model grid points. The 33x33 array o£ 


anomalous temperatures at the ten levels of the model for 
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September 1976 was then added to the model generated clima- 
tology for September 15 to constitute the initial tempera- 
meume conditions of the model; Similar 33x33x10 arrays for 
Wetober 1976 to Beary 1077 served as the verifying data 
PorecOmpucing verification statistics for the different 


model integrations as described below. 
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PEt. MODELING APPROACH 


In order to test a variety of atmospheric and oceanic 
processes that may be responsible for the formation and 
evolution of the large-scale temperature anomalies that are 
of interest, a series of initial value model integrations 
have been carried out using a 1lO-level primitive equation 
model in a closed rectangular basin. The eae LS that o£ 
Haney (1974) which has been improved to include time de- 
pendent seasonal forcing by the atmosphere, a parameteriza- 
tion of surface wind and convective mixing (Haney and 
Davies, 1976), and nonlinear lateral eddy viscosity based 
On two-dimensional turbulence theory (Haney and Wright, 
975). The ocean model is driven by wind and differential 
heating which is calculated from the model sea-surface 
temperature and prescribed values of the atmospheric solar 
radiation, total cloudiness, surface air temperature, rela- 
tive humidity and winds. These quantities are made up of 
a climatological part and in some cases an anomalous part. 

To initialize the model for a prediction experiment, 
the model is first integrated over a long "Spin-up" period 
mao) vyears) using climatological atmospheric forcing. The 
last year of this integration is referred to as the "model 
climatology" and the resulting currents and temperatures, 
especially in the mid-ocean areas, bear a close resemblance 


Gomtcnose in the North Pacific Ocean. Mies Lait lak -condLtewvons 
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for the prediction experiments are then obtained by combining 
a given anomaly, whether oceanic or atmospheric, with the 
model climatology for the same parameter for the appropriate 
time of the year. 

Once the anomalous atmospheric forcing and temperature 
anomaly fields were prepared for input into the model, five 
runs, each with different model configurations, were made. 
As stated in the introduction, the purpose of the different 
model experiments was to examine the relative importance of 
various atmospheric and oceanic processes in the formation 
and evolution of the observed anomalies. 

The five different model integrations are most easily de- 
Secribed by writing the equation of horizontal motion and the 


thermodynamic energy equation as follows, (Symbols defined in 
Table I), 


av aV 
He : le 
Sunk 3 oT 
pa tO, ®) +V°(AVT) + 5g (Kk FR) (2) 


The terms on the right hand side of (1) are the pressure 
gradient force, the Coriolis force, horizontal momentum dif- 
fusion and vertical diffusion. The terms on the right hand 
Side of (2) are the vertical mixing of heat due to the action 
Guesurrace wind and convection, horizontal diffusion and 
vertical diffusion. The ocean model is driven by a downward 
Flux of heat and momentum (stress) at the surface which 
emeers through boundary conditions in the vertical diffusion 


terms: 
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(K 99) GOwnwardsneat Elux at the surface 


N 
I 
© 
ll 
tO 
ll 


surface stress 


(K =>) = T 


In all of the model experiments climatological values of the 


heat flux, 0 ,were used. The stress, 1 ,was calculated from 


T= 70.C,|Vel Vo. 


where V< is a prescribed surface wind whieh 1s the sum of 
a climatological part and an anomalous part; the latter 
being calculated from Namias' sea level pressure anomalies 
described above. MThus, only Ve has an anomalous part, say 


ie and it enters the model equations in two places; as a 


memndary COndition (wind stress) in the vertical diffusion 


a 


mermem, (K a2) and directly in the vertical mixing term 


ea 
M(Q,T). The ability exists in the model then to allow the 
anomalous surface winds to enter into either, both or 
neither of the terms. In this way, we examine separately 
the role of anomalous wind driven currents and _ anomalous 
wind generated vertical mixing in the generation of oceanic 
thermal anomalies. 

In physical terms, the first model configuration, model 
1, describes the free evolution of an initial temperature 
anomaly since peo = 0 in this model. in model y2, ‘ehe 
effect of the anomalous wind stress in driving anomalous 


Ekman type flow is added, since Veo = Q only in the mixing 


berm in this model. According to steady Ekman theory this 
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Bow would be 9O" togthe right Of the anomalous surface 
winds. Model 3 incorporates the effects of anomalous wind 
stress in both Ekman type flow and surface wind mixing. 
This mixing process is expected to produce increased verti- 
cal mixing in regions of strong anomalous winds. Model 4 
is similar in its dynamics to model 2 but has the initial 
anomalous temperature fields equal to zero. Likewise, 
model 5 is similar to model 3 but with the initial anomalous 
temperatures equal to zero. These last two experiments are 
designed to evaluate the importance of having good initial 
conditions when making a prediction. 

In order to provide a quantitative comparison between 
the different models, a variety of statistics were computed 
for the five model cases at 30-day intervals of time. The 
Statistics calculated included the average absolute predic- 
tion error, root Mean square error and an Sl score associated 
with the model generated temperature anomalies. The follow- 
ing definitions indicate exactly what was measured by the 
statistics. Average absolute error and root mean square 
error are defined as follows: 1f ina series of N forecasts 
(lee. N different gridpoints), P represents the i-th fore- 
Gast and O. the corresponding observation, the average 


L 


absolute error is given by 
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and the root mean Square error by 


2 
¥ (F - 0.) 


RMSE = y 


(Brier and Alen, 1951). 

The Sl score is a dimensionless measure of the difference 
between predicted and observed horizontal gradients and is 
given by 
Zle.| 


ne 


(Teweles and Wobus, 1954), 
where re ls the error in the forecast temperature anomaly 
Seadient in the horizontal and G, is the observed tempera- 
ture anomaly gradient. The Sl score was further divided in 
order to measure the skill in forecasting both the north- 
south and the east-west gradient of the temperature anomaly. 
The various model predicted and TRANSPAC derived tem- 
perature anomaly fields were compared beginning with the 
October 1976 fields at 30-day intervals of time and at nine 
levels of the model. Statistics were computed only over an 
eae trom 35°N to 45°N and 130°W to 180°. The verification 
was restricted to the region east of the dateline because 
the region west of the dateline is generally characterized 
by numerous mesoscale temperature features (Bernstein and 
Meare, L97/) 2.5 .Sance thesintent of this investigation was in 
part to study present capabilities in modeling large scale 
thermal anomalies, it was assumed the model could not effec- 


tively predict such mesoscale features. 





IV. RESULTS AND CONCLUSIONS 


A vast amount of data is made available by this study 
and as is always the case, interpreting such a volume of 
watd im light of the original intentions of the study is 
Mis ticult. Describing the actual evolution of the thermal 
anomalies in the observed and various model cases is perhaps 
the most logical starting point. A later more detailed 
treatment of the verification statistics themselves will 
serve to highlight various aspects of the described results. 

Since the effects of anomalous wind stress in producing 
large scale thermal anomalies is considered paramount, a 
description of the anomalous atmospheric forcing for 
September 1976 to January 1977 should precede any descrip- 
tion of observed or model evolutions of thermal anomalies. 
Figure 2 shows the anomalous sea-level pressure fields for 
these months, contoured at 2 mb intervals with shaded areas 
indicating negative anomalous pressures. These anomalous 
pressure fields are shown only over the area of the tem- 
aoaaeure flelds given in figures 3 through 8. The September 
meomean cirreulation over the central Pacific was dominated 
myea trough at 700° mb which had progressed to a position 
south of the Aleutians. Downstream, the western North 
American ridge became established as an Arctic ridge 
(Taubensee, 1976). This mean 700 mb trough in the central 


Pacific is evidenced at the surface by negative anomalous 
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pressure. In October 1976, the mean 700 mb trough over the 
eastern Pacific retrograded somewhat, as did the downstream 
ridge over western North America (Wagner, 1977a). The nega- 
tive anomalous pressure pattern at sea level iS Tseen ern 
figure 2 to move westward in October, being replaced to the 
east by positive anomalous pressure associated with the 
ridge over western North America. The mean 700 mb trough 
Over the central Pacific remained relatively stationary but 
Mieensiflied Significantly (2 standard deviations from the 
mean) from November 1976 to January 1977 (Dickson, 1977; 
Taubensee, 1977; Wagner, 19775). This intensification is mani- 
fested in the sea level pressure by an extreme deepening 
of the negative anomalous pressure fields for November 
eacough January. 

The observed anomalous temperature data for September 
meee, November 1976 and January 1977 at 10, 102, 162 and 262 
meters are shown in figure 3. The observed evolution during 
this winter period in the upper levei is characterized by 
extreme cooling over all but the easternmost edge of the 
area. In the mid-levels (102 and 162 meters), the September 
field shows an extremely intense disturbance in the monthly 
Mean temperature anomaly at and just east of the dateline. 
This disturbance propagates westward out of the area by 
November. The 102 meter level exhibits a general cooling 
Over the winter months, Similar to that in the upper level. 
At 162 meters there is a general trend towards warming with 


the axis of temperature anomalies shifting from generally 


23 





north-south to east-west. The 262 meter field shows a less 

intense anomaly structure than the mid-levels for September 

but exhibits a more intense warming along the east-west axis 
than had occurred by January at 162 meters. 

Figure 4 shows the anomaly evolution predicted by model 1 
(climatological winds). The dominant feature in the upper 
levels is the development of an intense, basically mesoscale 
feature in the western portion of the area. This development 
is due to the strong north-south advection of mean tempera- 
ture by the anomalous geostrophic surface currents induced 
in the model by the intense disturbance in the western part 
of the basin at mid-levels. This same type of development 
occurs in the middle and lower levels so that a similar but 
less intense structure develops in the vertical. These 
features tend toward a slow westward propagation during the 
Winter months, indicating partial development and movement 
as free Rossby waves. 

Model 2 introduces anomalous winds, but only in the sur- 
face stress term (climatological winds are used in the verti- 
cal mixing term). Including this effect produces anomalous 
Ekman type flow in the upper levels which definitely improves 
the model predictions by cooling the upper layers (see figure 
5). The intense mesoscale features in the west produced by 
the strong north-south advection of mean temperature de- 
Scribed in model 1 are still very apparent in the anomaly 
Structure however. These anomalously cold mesoscale features 


are intensified and expanded in their horizontal scale in the 
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upper levels by the superposition of the large scale cooling 
resulting from the anomalous wind stress. In the middle to 
lower levels there iS a general warming, Similar to the ob- 
served case, which tends to reduce the intensity of the cold 
eddies. 

Model 3 includes the effects of anomalous winds in both 
the surface stress term and the vertical wind mixing of heat. 
This wind mixing process would be expected to produce in- 
creased vertical mixing in regions of increased surface winds 
pile as Seen in figure 6, the results for model 3 do not vary 
Significantly from model 2. In the ocean model, surface wind 
mixing depends in part on the net surface heat flux. If the 
net flux is downward, wind mixing occurs, but if the net sur- 
Face heat flux is upward as it is in this winter case there 
1s little or no wind mixing. Model 3 results then show only 
slight differences from model 2 in the uppermost layers. 

Model 4 includes anomalous winds only in the surface 
stress term (like model 2) but the initial temperature anomaly 
is set to zero in order to evaluate the importance of knowing 
eiewinitial conditions. In this particular winter case of 
very strong anomalous wind forcing, the temperature anomalies 
in the upper levels are predicted slightly better than in 
models 2 or 3, with general cooling over most of the area as 
Seen in figure 7. Results in middle to lower levels are 
worse, however, with the 102 m level showing development of 
small scale warm anomalies near 30-35N where cooling actually 


Pecurred and the 262 m level showing cooling to the northeast 
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where warming is observed. It may be noted at this point 
that the superposition of the results for model 4 with 
initial anomalous temperature equal zero and model 1 with 

no anomalous wind forcing will give almost exactly the 
results observed in model 2. This is cldéar evidence that 

in the model, the effects of external (wind) forcing and the 
effects of internal dynamical development take place rather 
independent of one another. 

Model 5 is the same in its dynamics as model 3 but again 
the initial anomalous temperature is set to zero. As pre- 
viously discussed, models 2 and 3 show little variation and 
meameseen in figure 8, model 5 results vary only slightly from 
model 4. 

The actual verification statistics, of course, provide more 
concrete evidence to the similarity and differences between 
the various models in successfully describing the evolution 
of the large scale temperature anomalies as observed. Table 
II presents these statistics in the following format: nine 
levels from 10 meters to 1700 meters depth are considered; 
the ratio of the root mean square error (RMSE) to the RMSE 
of climatology (a prediction of zero temperature anomaly 
everywhere) is given as the first number followed by a 
Similarly non-dimensionalized Sl score for 30 day intervals 
Eeom October 1976 to January 1977. A ratio of RMSE's and Sil's 
were used as a measure Of skill for easy interpretation since 
@ ratio less than one indicates skill over climatology 


whereas a ratio greater than one indicates less skill than 
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climatology. The skill of a forecast of persistence is also 
shown. These ratios for 10, 102, 162, and 262 meters have 
been plotted as a function of time in figures 9a and 9b for 
easy reference. 

The plot of skill as a function of time at 10 meters 
seen in figure 9a indicates an extremely rapid loss of skill 
as early as 30 days, followed by improved skill in the models 
Geeeo, 4 and 5) having wind forcing. Model 1 with no 
anomalous wind forcing shows no skill after 30 days but 
models 2 and 3 clearly demonstrate the benefits of consider- 
@mnG such anomalous forcing. Surprisingly, models 4 and 5 
with zero initial anomalous temperature conditions show con- 
stant improvement with time even over models 2 and 3. This 
would apparently indicate the greater importance of knowing 
the anomalous wind forcing than of knowing the initial tem- 
perature field exactly, especially in a case as this one 
where the anomalous winds were very strong. 

Perhaps the most conspicuous result is the fact that the 
model predictions at the surface are better when the initial 
anomaly is zero (model 4 and 5) than when it is the observed 
anomaly (model 2 and 3). In the latter case, strong geo- 
eemopiitc Ellow in the north=south direction 1s induced at 
the surface and this flow advects the mean isotherms in the 
model (cold to the north and warm *o the south) to produce 
a large amplitude disturbance in the western third of the 
basin which does not exist in the observed pattern. As de- 


Scribed earlier, the strong meridional flow is due to the 
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intense anomaly in the initial temperature field at and be- 
low 102 meters. One possible explanation for the model 
discrepancy is that the model simply does not propagate this 
developing disturbance westward and out of the domain fast 
enough. However, an examination of the surface anomaly 

data west of the dateline shows no large disturbance there 

to which the predicted disturbance in figure 4 could corre- 
Seong if the model error was simply one of phase error. 

A second possibility is that in September there existed a 
compensating deep temperature anomaly (below 400 meters) 
whose induced vertical shear nearly cancels that induced by 
the anomaly observed between 102 and 262 meters. in has 
case, the model assumption that the initial anomaly decreases 
Slowly to zero at depth, without changing sign (see DATA 
PROCESSING) is wrong. If this second explanation is true, 

1t implies that in order to predict or understand sea surface 
temperature anomalies, one must measure the anomaly field 

all the way down to the bottom of the sea. UntoOretunately,; 
imomwill require many case studies of this type before this 
question is adequately resolved. 

The plot of skill versus time for 102 meters in figure 
9a shows a somewhat less rapid loss of skill initially but 
demonstrates a fairly high amplitude oscillation about 1.0 
mieecene 30 to 120 day period of time. Comparison of the 
Curves for models 4 and 5 and models 1, 2 and 3 indicates 
a relation between these oscillations and the initial tem- 
perature conditions since this alternating behavior is not 


displayed by models 4 and 5. 





Figure 9b shows an extreme and rapid loss of skill (com- 
mBarsed to Climatology i.e. a forecast of zero anomaly) for 
persistence and the models at 162 and 262 meters. An oscil- 
lation in performance characterizes the 162 meter graph. 

It is important to note also that while models 4 and 5 have 
greater skill, they are not aS good as climatology. The 262 
meter plot does not exhibit this oscillatory characteristic, 
presumably as a result of less intense anomalous initial 
temperatures at this depth. It is rather disheartening that 
the model loses skill so fast at these levels. It 1s sur- 
prising too that persistence is also a very bad forecast 

at these depths which are generally believed to undergo 

Slow changes. 

The Sl scores computed to measure skill in predicting the 
gradient of large scale thermal anomalies are shown in Table 
I. Models 4 and 5 consistently demonstrate greater skill 
than either persistence or the other models in predicting 
gradients. This indicates of course that the ocean model has 
little skill in computing gradients since any field more 
closely resembling a zero field possesses greater skill. lis a 
the model Sl scores are broken up into north-south and east- 
west components, however, there is greater relative skill in 
Predicting north-south gradients than east-west gradients. 

| The results of this study can be summarized as follows. 
In the upper levels, the models which include anomalous wind 
forcing show skill over climatology out to 4 months, i.e. the 


ratio of RMSE of the model to RMSE of climatology is less 
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than one. However, Knowledge of the initial anomalous tem- 
perature conditions degrades the prediction skill. Two 
possible explanations were offered but a variety of addi- 
tional case studies are needed to properly resolve this 
problen. In the lower layers (162 to 262 meters), none of 
the models show skill over climatology and consideration 

of neither anomalous wind forcing nor initial anomalous tem- 
peratures improves the models predictive capabilities. The 
rather poor model results at these depths are a serious con- 
Semm because it is not likely that the addition of anomalous 
heat fluxes to the model predictions will have much effect 
at these levels. in -bhe siucure, 1t is planned to repeat 
this type of study for cases in which the anomalous atmo- 
spheric forcing was very weak and also to repeat the present 
study adding anomalous heat fluxes being computed by other 


NORPAX investigations. 
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V = NOrIizontal velocity 


H 
oa = reference density 
Pp = pressure 
fe = Coriolis parameter 
A = Worl ZOncal ditriusion Cockticient 
K = vertical diffusion coefficient 
dy = temperature 
M = parameterized vertical mixing term 
fe) = surface heat flux 
15 = surface wind stress 
TABLE I 
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